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1.0  INTRODUCTION 

After  an  icing  test  conducted  in  the  Arnold  Engineering  Development  Center  (AEDC) 
Engine  Test  Facility  in  1980,  a  study  was  initiated  to  obtain  a  better  understanding  of  the 
environmental  conditions  conducive  to  incipient  icing  on  the  leading  edge  of  an  airfoil.  The 
primary  result  was  the  development  of  a  steady-state  thermal  model  for  analyzing  the 
leading  edges  of  wings,  fins,  or  inlets  of  small  missiles  or  drones  at  subsonic  flight 
conditions.  The  model  was  found  to  be  helpful  in  predicting  the  effects  of  liquid  water 
content  (LWC)  on  the  leading-edge  surface  temperature  and  the  free-stream  static 
temperature  and  in  analyzing  incipient  icing  test  results. 

The  purpose  of  this  report  is  to  present  a  description  of  the  thermal  model,  the  model 
validation,  a  discussion  of  the  model  behavior,  and  examples  of  model  applications. 

2.0  MODEL  DEVELOPMENT 
2.1  APPROACH  AND  ASSUMPTIONS 

The  approach  was  to  design  a  steady-state  model  for  subsequent  trimming  and  validation 
based  on  experimental  or  test  results.  After  validation  of  the  steady-state  model,  further 
development  can  be  effected  to  improve  the  model  capability  for  addressing  a  wider  range  of 
conditions  and  time-dependent  input  parameters. 

Consistent  with  this  approach,  three  basic  assumptions  were  made  for  the  initial  model. 
First,  the  thermal  model  was  designed  for  steady-state  conditions,  i.e.,  when  thermal 
equilibrium  has  been  established.  In  flight  this  is  not  the  case  for  incipient  icing  during  the 
first  few  minutes  after  entry  into  an  icing  cloud,  when  the  airfoil  temperature  is  changing, 
but  it  is  the  case  after  the  temperature  has  stabilized.  Therefore,  one  needs  to  exercise  care  in 
interpreting  the  model  results  when  using  test  or  experimental  data  as  inputs.  Often  the  time- 
dependent  changes  in  test  conditions  are  sufficiently  small  to  permit  the  calculation  of  useful 
information.  Second,  the  model  does  not  account  for  ice  accretion.  Finally,  the  leading-edge 
shape  was  chosen  to  be  represented  by  a  half  cylinder  because  this  configuration 
approximates  typical  leading-edge  shapes  and  because  water  catch  efficiency  and  forced 
convection  heat-transfer  coefficient  data  for  a  cylinder  are  readily  available. 

Other  assumptions  and  considerations  will  be  discussed  in  this  section  and  subsequent 
sections  of  this  report. 
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2.2  DESCRIPTION 

The  thermal  model  is  a  heat  balance  of  the  leading  edge.  Based  on  the  conservation  of 
energy,  the  terms  of  the  heat  balance  equation  proffered  by  Messinger  (Ref.  1)  and  Schultz 
and  Wiilbanks  (Ref.  2)  were  used: 

QNET  -  Heat  Loss  -  Heat  Gain 

=  £|0SS  [QC  (Convective),  QE  (Evaporative),  QS  (Sensible)] 

-  Egain  [QF  (Air  Kinetic  Energy),  QKE  (Water  Kinetic  Energy)] 

To  this  heat  balance,  two  terms  were  added:  (1)  a  heat  loss  term  (QK,  Conduction)  to 
account  for  conduction  heat  loss  at  the  aft  structural  interface,  and  (2)  a  heat  gain  term  (PA, 
Heater  Power  Density)  to  account  for  an  anti-icing  heater.  These  heat-transfer  modes  are 
shown  schematically  in  Fig.  J .  The  basic  thermal  balance  equation  (per  unit  area)  is: 

QNET  =  HC(TS  -  TO)  +  MPR(M/A)L  +  (M/A)  (CW)  (TS  -  TO)  +  QK 
-  HC(R)V2/2GJ(CA)  -  (M/A)V2/2GJ  -  PA,  Btu/hr-ft2 

where 

A  =  Leading-edge  area,  ft2 

CA  =  Specific  heat  at  constant  pressure  for  air,  0.24  Btu/lbm-°F 
CW  =  Specific  heat  of  water,  1.0  Btu/lbm-°F 
G  =  Constant,  32.174  lbm-ft/lbf-sec2 
HC  =  Calculated  convective  heat-transfer  coefficient,  Btu/hr-ft2-°F 
J  =  Mechanical  equivalent  of  heat,  778  ft-lbf/Btu 
L  =  Latent  heat  of  vaporization  for  water,  1,060  Btu/lbm 
M  =  Rate  of  water  catch,  Ibm/hr 
MPR  =  Fraction  of  water  evaporated 
PA  =  Heater  power  density,  Btu/hr-ft2 
QK  =  Aft  conductivity,  Btu/hr-ft2 

R  =  Boundary-layer  recovery  factor,  0.875  average  for  laminar  and  turbulent 
flow  (Ref.  1).  (Use  the  appropriate  value(s)  based  on  the  local  flow 
conditions.) 

TO  =  Free-stream  static  temperature,  °F,  °R 
TS  =  Leading-edge  surface  (skin)  temperature,  °F,  °R 
V  =  Free-stream  velocity,  ft/sec 
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References  1  and  2  do  not  include  the  equations  for  calculating  HC  and  M,  (1)  and  (2) 
below.  The  following  equation  for  HC  was  obtained  from  Ref.  3: 

HC  =  0.194(TF)049  [V(RHO)/D]° 50  (1  -  |0/9O|3),  Btu/hr-ft2-°F  (1) 


where 

TF  -  0.5  (TS  +  TO),  °R 

RHO  =  Free-stream  air  density,  P0/R(TF),  lbm/ft3 
D  =  Effective  diameter  of  leading  edge,  ft  (See  Fig.  1) 

8  =  Included  angle  between  point  on  surface  and  stagnation  point,  integrated 
from  0  to  90  deg  to  obtain  an  average  value  over  the  leading  edge. 

(An  alternate  equation  for  HC,  from  Ref.  4,  was  used  for  verification: 

HC  =  (0.855)  (K/D)  (27,798.336  MU/K)0-4  [V(D)RHO/(MU)G]0-5,  Btu/hr-ft2-°F 

where 

K  =  Conduction  heat-transfer  coefficient  for  air  in  the  boundary  layer 
=  0.001533  TFK'  V(TFK  +  245.4  x  10-12/TFK),  Btu/hr-ft-°F  (Ref.  5) 

(TFK  =  TF,  °K), 

and 

MU  =  Viscosity  of  air  in  boundary  layer 
=  lbf-sec/ft2,  calculated  using 
Sutherland’s  equation. 

Agreement  within  ±  1.1  percent  was  obtained.) 

M  =  0.38(VK)D(LWC)E(S),  lbm/hr  (2) 


where 

VK  =  Velocity,  kts 
LWC  =  Liquid  water  content,  gm/m3 
E  =  Water  catch  efficiency. 
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Excluding  the  S  term,  the  aforementioned  equation  for  M  as  obtained  from  Refs.  4  and  6 
gives  the  water  catch  rate  per  foot  of  airfoil  span.  The  S  term  was  added  to  account  for  the 
total  span  in  feet. 

The  water  catch  efficiency,  E,  is  the  ratio  of  the  water  intercepted  (impinged  on  the 
leading  edge)  to  the  amount  of  water  contained  in  the  cloud  streamtube  which  flows  over  the 
leading  edge.  It  is  a  function  of  the  modified  inertia  parameter,  KO,  which  accounts  for 
water  droplet  inertia.  KO  is  calculated  from  the  expression  given  in  Ref.  4  (accurate  within  ± 
5  percent,  per  Ref.  6): 

KO  =  1.87  x  10-7  [1.15  VK/(MU)G]0-6  DD1  6/12  D(RHO)0-4 


where 

DD  =  Volume  median  droplet  diam,  p. 

Water  catch  efficiency  is  determined  from  a  graph  of  E  as  a  function  of  KO  for  a  cylinder  as 
obtained  from  Ref.  6.  This  relationship  is  shown  in  Fig.  2  in  addition  to  that  for  other 
aerodynamic  configurations. 

The  next  term  in  the  thermal  balance  to  be  defined  is  the  evaporation  fraction,  MPR,  the 
ratio  of  the  amount  of  water  evaporated  from  the  surface  to  the  amount  of  water  intercepted 
by  the  surface,  or 


MPR  -  EP/ (M/A) 


where 

EP  =  Evaporation  potential,  for  this  analysis  the  potential  or  maximum  amount 
of  water  per  unit  area  that  can  be  evaporated,  a  function  of  the  difference 
between  the  water  vapor  pressure  based  on  the  surface  (leading  edge) 
temperature  and  the  vapor  pressure  based  on  the  free-stream  static 
temperature. 

The  intercepted  water  is  subject  to  three  processes:  freezing,  being  windswept  along  the 
surface  (runback),  and  evaporation.  In  this  analysis,  no  provisions  are  made  to  account  for 
freezing  or  ice  accretion.  Therefore,  only  complete  evaporation  or  partial  evaporation  and 
partial  runback  can  occur.  Complete  evaporation  can  be  effected  only  at  relatively  high 
surface  temperatures,  generally  above  100°F  for  the  usual  cloud  encounters.  Therefore,  the 
more  common  condition  is  partial  evaporation  of  the  intercepted  water  with  attendant 
runback  of  the  remainder.  However,  the  partial  evaporation  in  this  model  is  the  maximum 
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thal  can  be  potentially  evaporated  based  on  the  two  vapor  pressures  previously  defined.  This 
is  the  general  approach  used  by  Messinger,  and 

EP  =  KM  (HC)  (PTSW  -  PT0W)/P0,  lbm/hr-ft2 


where 

PO  =  Free-stream  static  (ambient  or  atmospheric)  pressure,  psia 

PTSW  =  Water  vapor  pressure  at  the  surface  temperature,  TS.  For  TS  >  32°F, 
PTSW  was  obtained  from  a  curve  fit  of  vapor  pressure  as  a  function  of 
temperature  from  steam  tables  (Ref.  7). 

PTOW  =  Water  vapor  pressure  at  TO,  i.e. ,  free-stream  (atmospheric)  vapor  pressure 
of  the  water  droplets,  psia.  This  assumes  that  the  water  droplet 
temperature  is  the  same  as  the  free-stream  static  temperature.  For 
TO  <  32°F,  the  following  expression  was  used  for  PTOW  (Refs.  2  and  8): 

PTOW'  -  14.696  exp  (2.3  [A1  +  Bl/T  +  C1(T2-K)/T  [lODi(T2-K)2_  i<0] 

+  El  [10^4.11-1^%  pSia 

where 

T  =  (TO  4-  459. 67)/ 1.8,  °K 
and 


A1  = 

5.4266514 

D1  = 

1.1965  x  10-" 

K  =  2.937  x  105 

B1  = 

-2005.1 

El  - 

-4.4  x  10~3 

Cl  = 

1.3869  x  104 

FI  = 

-5.7148  x  10-3 

For  TO  >32°F,  PTOW'  was  obtained  from  a  curve  fit  of  vapor  pressure  as  a  function  of 
temperature  from  steam  tables  (Ref.  7). 

Referring  to  the  EP  equation,  KM  according  to  Messinger  is  a  constant,  2.9.  Sogin  (Ref. 
9)  derives  it  from  the  expression  0.622 1/CA  =  2,90,  i.e.,  where  I  =  (PR/SQ0-67  =1.12  and 
states  that  it  is  based  on  the  assumption  of  laminar  mass  transfer,  the  mass  transfer  from  the 
surface  to  the  free  stream.  Sogin  further  states  that  KM  =  2.6  if  the  mass  transfer  is 
turbulent  (I  =  1 ,0).  He  includes  other  estimates  of  I  by  many  investigators  for  the  reader’s 
perusal.  In  this  model  KM  is  calculated  for  each  environmental  condition  using  the 
expression 
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KM  =  (0.622/C A)  (PR/SC)0-67,  lbm-°F/Btu 


where 


PR  =  3600(G)CA(MU)/K  =  27798.336  MU/K 


and 


SC  =  MU/RHO(DF) 


where 


DF  =  0.000146  TF2  5/((TF  +  441)(PO/14.696)J,  ftVhr  (Ref.  10) 

KM  was  evaluated  for  a  set  of  72  data  points  and  varied  from  2.8602  to  2.8686  (I  =  1 .104  to 
1.107). 

The  ratio  EP/(M/A)  or  MPR  is  an  indicator  of  the  amount  of  runback.  When  EP 
=  M/A, 

MPR  =  1.0  (no  runback,  all  intercepted  water  is  evaporated). 

When  EP  <  M/A, 

MPR  =  EP/(M/A).  (M/A  -  EP  =  unit  runback). 

Referring  again  to  the  EP  equation,  some  authors  include  an  additional  term  called  the 
“wettedness”  factor,  particularly  for  analyzing  the  wing  areas  aft  of  the  leading  edge.  This 
factor  is  the  ratio  of  the  wetted  surface  area  to  the  total  area.  For  the  leading  edge,  the 
wettedness  factor  is  generally  assumed  to  be  unity.  Such  is  the  case  in  this  model.  The  reader 
is  cautioned  to  use  a  value  appropriate  to  the  particular  conditions/configuration  being 
investigated. 

The  heat  gain  (QG)  and  the  heat  loss  (QL)  are  calculated  as  follows: 

QG  =  QF  +  QKE  +  PA,  w/in.2 
QL  =  QC  +  QE  +  QS  +  QK,  w/in.2 
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where  all  of  the  component  terms  are  converted  from  Btu/hr-ft2  to  w/in.2  by  the  factor 
2.0353  x  IO-3. 


QNET  =  QG  -  QL,  w/in.2 

The  units  w/in.2  are  used  because  they  are  the  commonly  used  units  for  describing  the  heater 
capacity  or  the  heater  power  density. 

3.0  MODEL  VALIDATION 


3.1  APPROACH 


The  model  validation  consisted  of  inputting  measured  test  data  to  the  thermal  model, 
examining  the  degree  of  thermal  balance  achieved,  and  then  adjusting  the  model  to  optimize 
the  thermal  balance.  For  an  ideal  thermal  balance,  the  heat  gain  (QG)  should  equal  the  heat 
loss  (QL).  The  mode!  requires  the  following  inputs: 


Free-stream  static  pressure  (altitude  ambient),  P0 

Free-stream  static  temperature,  TO 

Flight  Mach  number,  XMO 

Leading-edge  surface  temperature,  TSM 

Leading-edge  surface  area,  A 

Airfoil  span  (inlet  perimeter),  S 

Leading-edge  thickness,  D 

Free-stream  liquid  water  content  (cloud),  LWC 

Volume  median  water  droplet  diameter,  DD 

Leading-edge  heater  power  density,  PA 

The  model  outputs  are: 

Recovery  of  the  water  kinetic  energy,  QKE 
Recovery  of  the  air  kinetic  energy,  QF 
Total  heat  gain,  QG  (QKE  +  QF  +  PA) 

Convective  heat  transfer,  QC 
Evaporative  heat  transfer,  QE 
Sensible  heat  transfer,  QS 

Conduction  heat  transfer,  QK  (between  the  leading  edge  and  aft  structural  interface) 
Total  heat  loss,  QL  (QC  +  QE  +  QS  +  QK) 

Net  heat  transfer,  QNET  (QG  —  QL) 

Calculated  surface  temperature,  TSC 
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The  calculated  data  should  result  in  thermal  balance  for  subsonic  flight  conditions 
provided  that  the  input  data  are  correct,  the  model  is  accurate,  and  the  leading  edge  has 
reached  thermal  equilibrium  (steady-state  conditions).  Two  indicators  were  used  to  evaluate 
the  degree  of  thermal  balance  achieved.  The  first  indicator  was  the  dependent  parameter 
QNET,  which  should  equal  zero;  i.e.,  the  sum  of  the  heat  gain  terms  (QG  =  QKE  +  QF  + 
PA)  should  equal  the  sum  of  the  heat  loss  terms  (QL  =  QC  +  QE  +  QS  +  QK),  or  QNET 
=  QG  -  QL  -  0.0.  The  second  was  the  leading-edge  surface  temperature;  i.e.,  the 
calculated  surface  temperature  (TSC)  should  equal  the  measured  surface  temperature  (TSM) 
when  thermal  balance  is  achieved.  The  calculated  surface  temperature  was  obtained  by 
iteration  of  the  model  on  the  surface  temperature  until  QNET  *  0.0.  (In  the  model, 
iteration  is  complete  when  |  QNET  |  <  0.01  w/in.2).  Since  an  exact  thermal  balance  was  not 
achieved  on  a  data-point-by-data-point  basis,  the  magnitude  of  the  difference  between  the 
mean  of  the  data  set  (linear  curve-fit  of  the  calculated  versus  measured  temperatures)  and 
the  line  of  ideal  surface  temperature  agreement  can  be  used  as  an  indicator  of  the  extent  of 
the  thermal  unbalance  in  terms  of  surface  temperature,  °F.  The  temperature  difference  in  qF 
is  a  more  comprehendible  gauge  of  the  degree  of  thermal  balance  achieved  than  is  the 
unbalance  in  QNET  in  w/in.2. 

The  available  data  used  for  the  validation  were  obtained  from  AEDC  tests  of  two  similar 
leading-edge  configurations,  A  and  B.  Both  configurations  are  characterized  by  the  leading- 
edge  schematic  shown  in  Fig.  1  where  the  thickness,  D,  was  0.22  and  0.54  in.  for 
configurations  A  and  B,  respectively.  Both  configurations  were  equipped  with  electrical 
resistance  anti-icing  heaters.  Only  Configuration  B  was  equipped  with  thermocouples  for 
measuring  the  temperature  gradient  across  the  aft  structural  interface.  The  gradients  were 
used  to  estimate  the  conduction  heat  loss. 

Both  dry-  and  wet-air  data  were  acquired  during  the  tests  of  both  configurations,  with 
and  without  the  heaters  in  operation.  The  simulation  of  cloud  icing  conditions  was 
accomplished  with  the  icing  water  spray  system.  The  water  spray  system  was  operated  to 
produce  droplets  with  a  mass  median  diameter  of  approximately  20  n  and  a  uniform 
distribution  of  droplets  throughout  the  free-jet  airflow  based  on  a  prior  calibration  of  the 
spray  system. 

The  usual  test  procedure  was  to  acquire  steady-state  dry-air  test  data  at  a  constant  upper 
airflow  temperature  with  the  heater  on  or  off  (as  required  to  meet  the  test  objective), 
activate  the  icing  water  spray  system  to  produce  inlet  icing  conditions  (data  system  off),  and 
then  gradually  decrease  the  airflow  temperature  1  to  2°F/min  while  recording  data  in  the 
transient  mode.  This  very  slow  change  in  conditions  was  considered  to  approximate  steady- 
state  data.  During  one  test,  however,  the  data  were  acquired  at  steady-state  conditions. 
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The  range  of  flight  conditions  of  the  validation  data  ranged  from  sea  level  to  13,300-ft 
altitude,  0.45  to  0.80  flight  Mach  No.,  -  6.0  to  35°F  free-stream  static  temperature,  and  0.2 
to  1.5  gm/m3  liquid  water  content.  Leading-edge  surface  temperatures  ranged  from  10  to 
80°F. 

The  validation  was  accomplished  using  two  groups  of  data  for  both  configurations,  dry- 
air  (cloud  off)  data  and  wet-air  (cloud  simulation)  data.  Because  of  data  scatter,  each  group 
is  treated  as  a  data  set,  and,  where  appropriate,  a  first-order  curve  fit  of  the  data  set  is 
compared  to  the  line  of  ideal  agreement.  The  estimate  of  the  one-sigma  standard  deviation 
(SD)  of  the  data  from  the  curve  fit  is  also  shown. 

In  the  first  group,  the  effects  of  conduction  heat  loss  across  the  aft  structural  interface 
are  discussed  including  the  procedure  used  to  compute  the  conduction  heat  loss.  The 
conduction  heat  loss  relationships  developed  during  the  dry-air  analysis  are  then  used  foT  the 
wet-air  data  analysis. 

3.2  THERMAL  BALANCE,  CONFIGURATION  A 
3.2.1  Thermal  Balance,  l>r\  Air  (Cloud  Off) 

For  the  dry-air  case,  the  wrater-dependent  terms  (QE,  QS,  QKE)  in  the  thermal  balance 
equation  (Fig.  1)  are  zero,  and  the  equation  reduces  to 

QNET  =  QC  +  QK  -  QF  -  PA 
Initially,  QK  was  assumed  to  be  zero.  Therefore,  the  equation  reduces  to 

QNET  =  QC  -  QF  -  PA 

The  dry-air  thermal  balance  achieved  for  the  Configuration  A  test  data  are  presented  in 
Fig.  3,  heater  on  and  off.  Ideally,  if  perfect  thermal  balance  were  achieved,  all  of  the  data 
would  lie  on  the  dashed  line.  Experimentally,  data  scatter  occurred;  therefore,  a  first-order 
curve  fit  was  drawn  through  the  data  set,  as  shown.  The  maximum  difference  between  the 
ideal  and  experimental  lines  was  about  0.3  w'/in.2  near  the  20-w/in.2  level  of  heat  loss  (or 
gain).  This  is  considered  good  agreement.  Since  thermal  balance  was  essentially  achieved, 
the  initial  assumption  that  QK  -  0  was  valid.  The  estimate  of  the  one-sigma  standard 
deviation  of  the  data  set  (SD)  from  the  curve-fit  line  was  ±1.23  w/in.2.  The  maximum 
deviation  for  a  single  point  from  the  ideal  agreement  line  w'as  2.0  w/in.2. 
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The  surface  temperature  data  for  the  same  data  set  are  shown  in  Fig.  4.  For  the  same 
measured  temperature,  the  calculated  surface  temperature  was  up  to  7.5°F  less  than  the 
temperature  required  for  thermal  balance  (corresponds  to  the  2  w/in.2  in  Fig.  3).  Again,  a 
curve  fit  for  the  data  set  was  drawn  and  differed  from  the  QNET  =  0.0  line  (ideal 
agreement)  by  -3°F  for  measured  surface  temperatures  near  50°F  (corresponds  to  0.3 
w/in.2.  Fig.  3)  and  +l.3°F  at  a  measured  surface  temperature  of  10°F.  This  good 
agreement  is  interpreted  to  mean  that  the  QC  and  QF  terms  in  the  thermal  model  are  valid. 

The  difference  between  the  measured  and  calculated  leading-edge  surface  temperatures  is 
presented  in  Fig.  5.  The  data  spread  was  about  13°F.  The  calculated  temperatures  were 
within  ±4.0°F  (SD)  of  the  measured  temperatures  for  the  data  set. 

3.2.2  Thermal  Balance,  Wet  Air  (Cloud  Simulation) 

For  the  wet-air  case,  all  of  the  terms  in  the  thermal  balance  equation  are  used  except  QK; 
again,  QK  was  assumed  to  be  zero.  Therefore,  the  equation  is 

QNET  =  QC  +  QE  +  QS  -  QF  -  QKE  -  PA 

The  thermal  balance  for  the  wet-air  data  set  (only  heater-on  data  were  available)  are 
presented  in  Fig.  6.  The  data  scatter  was  worse  than  the  dry-air  data  set,  and  the  curve  fit 
differed  from  the  ideal  agreement  line  by  7.7  w/in.2  at  the  heat  loss  (QL)  level  of  about  27 
w/in.2  The  difference  between  the  calculated  and  measured  surface  temperatures 
corresponding  to  the  same  data  set  for  the  same  27-w/in.2  heat  loss  level  was  about  9.0°  F 
(Fig.  7).  The  data  in  Fig.  6  indicate  that  the  heat  loss  was  too  high  (or  the  heat  gain  too  low) 
and  that  an  adjustment  of  the  heat  balance  equation  is  required. 

To  adjust  the  model,  the  terms  in  the  heat  balance  equation  were  evaluated  to  determine 
which  term  or  terms  might  be  adjusted  to  effect  thermal  balance  of  the  Configuration  A  wet- 
air  data  set.  The  evaporative  heat-transfer  term,  QE,  was  selected  as  the  term  for  adjustment 
based  on  the  following  rationale.  First  the  other  terms  in  the  equation  are  believed  to  be 
reasonably  accurate  because: 

1 .  QC  and  QF  were  validated  by  the  results  obtained  using  the  dry-air  data, 

2.  QS  and  QKE  are  relatively  straightforward  calculations  (Section  2.0),  and 

3.  PA  is  a  directly  measured  input. 
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Of  these  five  terms,  QF,  QKE,  and  PA  comprise  the  net  heat  gain;  therefore,  the  net  heat 
gain  must  be  approximately  correct.  Improvement  in  the  thermal  balance,  then,  can  be 
achieved  only  by  a  reduction  of  the  net  heat  loss  (Fig.  6).  Of  the  candidate  heat  loss  terms 
QC,  QS,  and  QE,  QC  and  QS  are  believed  to  be  reasonably  accurate  for  the  reasons  just 
given.  This  Leaves  QE  as  the  choice  for  adjustment  to  improve  the  thermal  balance. 
Furthermore,  the  mass  and  heat-transfer  mechanisms  related  to  the  evaporative  process  on 
the  leading  edge  are  believed  to  be  the  most  subject  to  error.  Continued  model  development 
may  reveal  a  more  exacting  adjustment  procedure. 

Continuing  with  the  QE  adjustment  discussion,  reduction  of  the  net  heat  loss  to  achieve 
thermal  balance  will  require  a  reduction  in  QE.  This  can  be  accomplished  by  multiplying  QE 
by  a  fractional  correlation  coefficient.  Application  of  a  correlation  coefficient  of  0.3  (0.3 
x  QE)  resulted  in  the  optimum  balance,  i.e. ,  positive  coefficients  greater  than  or  less  than 
0.3  degraded  the  thermal  balance.  The  results  are  shown  in  Fig.  8  where  the  difference 
between  the  ideal  line  and  the  data  set  curve-fit  line  is  about  1.8  w/in.2  compared  to  the 
unadjusted  difference  of  about  7.7  w/in.2  (Fig.  6).  The  unadjusted  maximum  temperature 
difference  was  23°F  (Fig.  7)  compared  to  the  adjusted  maximum  of  6°F  shown  in  Fig.  9. 
The  degree  of  thermal  balance  was  less  than  that  obtained  for  the  dry-air  data  set,  probably 
because  of  bias  in  the  model.  At  32°F  surface  temperature,  the  model  bias  was  +4.5°F. 

The  difference  between  the  measured  and  the  calculated  leading-edge  surface  temperatures 
for  the  Configuration  A  wet-air  data  are  shown  in  Fig.  10.  The  data  spread  and  the  standard 
deviation  were  greater  than  that  obtained  with  the  dry-air  data  set,  18°F  (compared  to  13°F, 
Fig.  5)  and  ±4.8°F  (compared  to  ±4°F,  Fig.  5).  Some  increase  in  the  scatter  was  expected 
because  the  water-dependent  terms  (QE,  QS,  and  QKE)  are  used  for  the  wet-air  data  but  are 
zero  for  the  dry-air  case. 

3.2,3  Data  Summary,  Configuration  A 

For  the  Configuration  A  dry-air  test  data  set,  the  thermal  model  predicted  the  surface 
temperatures  within  about  +1  to  —  3°F,  and,  for  the  wet-air  test  data  set,  within  about  +  5 
to  -6°F,  the  latter  after  adjustment  of  the  evaporative  heat-transfer  term,  QE.  The 
standard  deviations  of  the  differences  between  the  calculated  and  the  measured  leading-edge 
surface  temperatures  were  ±4  and  ±4.8°F  for  the  dry-air  and  the  wet-air  data  sets, 
respectively.  At  32°F,  the  model  bias  was  +  4.5°F  for  the  wet-air  data  set. 

Although  the  adjustment  to  QE  alone  did  improve  the  data  agreement,  continued 
analysis  and  experimentation  should  result  in  a  better  understanding  of  all  of  the  thermal 
and  mass  transfer  processes  and  an  improved  model. 
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3.3  THERMAL  BALANCE,  CONFIGURATION  B 

3.3.1  Thermal  Balance,  Dry  Air  (Cloud  Off) 

For  the  Configuration  B  dry-air  case,  the  water-dependent  terms  in  the  energy  balance 
equation  are  zero,  as  before.  The  initial  step  in  the  analysis  was  to  assume  that  QK  =  0.0 
and  to  plot  the  heat  balance  results  using  the  Configuration  B  dry-air  data  set.  These  dry-air 
data  were  the  only  available  data  acquired  at  steady-state  conditions.  The  results  in  Fig.  11 
show  an  unbalanced  condition  for  the  heater-on  data.  Since  the  heat  gain  was  too  high  by 
7.5  w/in.2,  the  unbalance  was  attributed  to  not  accounting  for  the  conduction  heat  loss.  The 
conduction  heat  loss  should  be  a  function  of  the  temperature  differential  across  the  aft 
structural  interface.  The  available  data  were  plotted  as  shown  in  Fig.  12,  and  a  linear  curve 
fit  was  derived  from  the  data  and  used  to  compute  the  conduction  heat  loss,  QK. 

The  model  was  rerun  with  the  conduction  heat  loss  calculation,  and  the  data  were 
replotted;  the  results  are  presented  in  Fig.  13.  The  data  were  within  2.0  w/in.2  of  the  ideal 
agreement  line,  and  the  curve  fit  of  the  data  was  coincident  with  the  ideal  agreement  line. 
The  corresponding  calculated  versus  measured  surface  temperature  data  set  is  shown  in  Fig. 
14;  again,  the  data  set  curve-fit  line  was  coincident  with  the  ideal  agreement  line.  The 
standard  deviation  for  the  data  set  curve  fit  was  ±  1.03°F. 

The  differences  between  the  measured  and  calculated  surface  temperature  data  for  the 
dry-air  conditions  are  shown  in  Fig.  15.  This  set  of  Configuration  B  data  exhibited  much 
better  agreement  than  did  the  Configuration  A  data  set.  Of  all' the  data  presented  in  this 
report,  this  data  set  was  the  only  data  acquired  under  steady-state  conditions.  The  standard 
deviation  for  this  steady-state  set  was  ±  1°F  compared  to  the  Configuration  A  dry-air  data 
standard  deviation  of  ±4QF. 

3.3.2  Thermal  Balance,  Wet  Air  (Cloud  Simulation) 

The  same  QK  calculation  derived  from  the  dry-air  data  was  also  used  for  the  wet-air  data 
analysis.  The  wet-air  data  computed  with  this  QK  calculation  are  shown  in  Fig.  16.  The  total 
heat  loss  exceeded  the  total  heat  gain;  therefore,  thermal  balance  was  not  achieved.  The 
maximum  difference  between  the  data  curve-fit  and  the  ideal  agreement  line  was  10  w/in.2. 
The  corresponding  surface  temperature  curve-fit  for  the  same  data  set  indicated  a  12° F 
difference  at  a  measured  surface  temperature  of  60°F  (Fig.  17). 

Since  a  similar  unbalance  was  experienced  with  the  Configuration  A  wet-air  data  set,  the 
Configuration  B  data  were  recomputed  using  the  same  0.3  correlation  coefficient  for  the 
evaporative  heat  loss  as  was  used  for  the  Configuration  A  data.  The  10  w/in.2  unbalance 
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(Fig.  16)  was  reduced  to  2.8  w/in.2  as  shown  in  Fig.  18,  a  much  improved  degree  of  thermal 
balance.  The  corresponding  surface  temperature  agreement  was  improved  from  a  maximum 
of  -  12°F  (Fig.  17)  to  that  shown  in  Fig.  19,  -3.5°F.  However,  the  standard  deviation  of 
the  data  was  high  (±3.04°F).  The  model  bias  at  32°F  measured  surface  temperature  was 
+  3.2°F. 

The  differences  between  the  measured  and  calculated  surface  temperature  data  for  the 
wet-air  conditions  are  presented  in  Fig.  20.  These  Configuration  B  data  exhibited  a  similar 
range  of  temperature  differences  with  about  a  ±  1°F  improvement  in  the  standard  deviation 
compared  to  the  Configuration  A  wet-air  data  set  (Fig.  10). 

3.3.3  Data  Summary,  Configuration  B 

After  accounting  for  aft  structural  interface  conduction  heat  loss,  excellent  agreement 
was  obtained  between  the  calculated  and  measured  surface  temperatures  for  the  dry-air  data 
set.  For  the  wet-air  data  set  computed  with  the  model  adjustment,  the  calculated-to- 
measured  surface  temperature  curve-fit  line  was  a  maximum  of  3.5°F  below  the  ideal 
agreement  line,  but  the  standard  deviation  of  the  set  was  high  (3.04°F).  The  standard 
deviations  of  the  differences  between  the  calculated  and  measured  leading-edge  surface 
temperatures  were  ±  1  and  ±4°F  for  the  dry-air  and  wet-air  data  sets,  respectively.  For  the 
wet-air  data  set,  the  model  bias  was  +3.2CF  at  32°F  measured  surface  temperature. 

3.4  DATA  UNCERTAINTY 

The  uncertainties  of  the  measured  test  data  used  as  input  parameters  to  the  model  were 
propagated  through  the  model  to  obtain  an  estimated  uncertainty  of  the  primary  calculated 
parameters,  heat  loss,  QL,  and  heat  gain,  QG.  The  sensitivities  (influence  coefficients)  for 
the  model  discussed  in  Section  4.1  of  this  report  were  combined  with  the  bias  and  precision 
of  the  measured  data  using  the  procedure  described  in  Ref.  1 1 .  The  results  for  a  nominal 
flight  condition  are  shown  in  Table  1.  The  uncertainty  of  the  model  is  not  know'n  but  is 
probably  no  better  than  the  5-  to  9-percent  uncertainties  shown  in  Table  1.  Therefore,  the 
total  uncertainty  of  the  calculated  data  obtained  from  the  measured  inputs  and  the  model  is 
probably  on  the  order  of  15  to  20  percent. 

3.5  SUMMARY 

Using  the  available  AEDC  test  data,  the  degree  of  thermal  balance  achieved  and  the 
corresponding  correlation  of  the  calculated  and  measured  leading-edge  surface  temperatures 
were  used  to  validate  the  thermal  model.  An  evaporative  heat  loss  correlation  coefficient, 
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EC,  of  0.3  was  used  to  optimize  the  agreement  between  the  calculated  and  measured  surface 
temperatures.  The  resultant  heat  balance  equation  thus  becomes 

QNET  =  HC(TS  -  TO)  +  (EC)  (MPR)  (M/A)L  +  (M/A)CW(TS  -  TO) 

+  QK  -  HC(R)VV2GJ(CA)  -  (M/A)W2GJ  -  PA 


where 


EC  =  0.3 

Using  this  coefficient,  the  model  is  considered  validated  to  the  extent  possible  with  the 
two  groups  of  available  data.  At  32°F  surface  temperature,  the  average  model  bias  is  about 
4°F.  Continued  thermal  model  development  should  improve  the  model  accuracy. 

4.0  MODEL  OUTPUT  PARAMETER  CHARACTERISTICS 

The  validated  model  (with  the  0.3  evaporative  correlation  coefficient)  was  operated  at 
typical  incipient  icing  conditions  with  the  input  parameters  varied  to  determine  the 
corresponding  effects  on  the  output  or  calculated  parameters.  The  results  are  presented  to 
assist  the  reader  in  evaluating  the  relative  effects.  The  aft  interface  conduction  heat  loss 
term,  QK,  was  omitted  to  simulate  an  adiabatic  leading-edge  structural  interface  (Fig.  1). 

4.1  MODEL  SENSITIVITY 

Each  input  parameter  was  varied  4- 1  percent  from  the  typical  flight  condition  of  Mach 
No.  0.55,  3,000-ft  altitude,  w  here  the  altitude  ambient  static  pressure,  PO,  was  based  on  Ref. 
12.  The  altitude  ambient  static  temperature  was  chosen  such  that  thermal  balance  is  achieved 
(i.e.,  QNET  =  0.0)  and  is  a  temperature  conducive  to  icing.  All  of  the  input  conditions 
chosen  for  this  typical  flight  condition  are  shown  in  the  third  vertical  column  in  Table  2  and 
are  designated  the  base  conditions.  Without  the  1 -percent  perturbation,  the  values  of  the 
base  output  parameters  are  shown  for  reference  on  the  second  line  of  the  table.  The 
remainder  of  the  table  lists  the  percent  changes  in  each  output  parameter  for  a  1 -percent 
increase  in  each  input  parameter.  The  +  1-percent  change  data  are  presented  in  Fig.  21  to 
show'  the  relative  effects  on  the  output  parameters,  which  generally  varied  from  - 1 .0  to 
+  1.0  percent.  However,  a  1-percent  increase  in  TS  resulted  in  a  1.86-percent  increase  in  the 
evaporative  heat  transfer  (QE).  A  1-percent  increase  in  XM0  resulted  in  a  2.5-percent 
increase  in  QF  and  a  3.1-percent  increase  in  QKE. 
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The  changes  in  QNET  in  w/in.2  are  shown  in  Fig.  22.  QNET  is  sensitive  to  XMO  but  is 
most  sensitive  to  TS.  Therefore,  when  preparing  the  criteria  for  a  test  or  experiment, 
accurate  measurement  of  these  parameters  is  essential. 

4.2  HEAT  LOSS  AUDIT 

The  relative  contribution  of  each  of  the  heat  loss  terms  (QC,  QE,  QS)  to  the  total  heat 
loss  (QL)  was  calculated  as  the  percent  of  the  total  heat  loss  and  is  presented  in  Fig.  23  as  a 
function  of  liquid  water  content  (LWC)  for  one  flight  condition.  At  zero  LWC,  all  of  the 
heat  loss  is  convective  (100  percent).  As  LWC  is  increased,  QC  decreases  and  QE  increases 
up  to  LWC  =  0.0165  gm/m3.  For  these  lower  portions  of  the  QC  and  QE  curves  (LWC  < 
0.0165),  all  of  the  intercepted  water  is  being  evaporated,  i.e.,  MPR  =  1.0.  Above  LWC 
=  0.0165  as  LWC  is  further  increased,  the  percentage  of  water  evaporated  progressively 
decreases  (MPR  <  1.0)— the  remainder  being  runback.  The  convective  heat  transfer,  QC, 
and  the  evaporative  heat  transfer,  QE,  decrease  as  the  percentage  of  sensible  heat  transfer, 
QS,  progressively  increases  with  increasing  LWC. 

The  calculated  QC  and  QE  data  exhibited  approximately  16  percentage  point  changes 
between  LWC  values  of  zero  and  0.0165  gm/m3.  This  vividly  indicates  that  the  maximum 
evaporative  cooling,  percentage-wise,  occurs  at  low  values  of  LWC. 

5.0  MODEL  APPLICATIONS 

The  model  may  be  used  to  analyze  the  thermal  energy  balance  on  the  leading  edge  of 
small  airfoils.  To  date,  application  of  this  model  has  been  limited  to  leading-edge  thickness 
(D)  values  near  0.5  in.,  volume  (or  mass)  median  droplet  diameters  (DD)  from  10  to  30 
microns,  altitudes  below  15,000  ft,  flight  Mach  numbers  (XMO)  up  to  0.8,  free-stream  static 
temperatures  (TO)  from  -25  to  32°F,  leading-edge  surface  temperatures  (TS)  from  32°F 
(instantaneously)  to  130°F  (continuously),  and  liquid  water  contents  (LWC)  up  to  3.0 
gm/m3. 

The  model  uses  an  estimate  of  the  boundary-layer  temperature  (TF)  as  suggested  in  Refs. 
4  and  13  which  is  the  average  of  the  free-stream  static  temperature  and  the  surface 
temperature  for  the  calculation  of  transport  properties  (e.g.,  HC)  and  the  associated  air 
properties  (MU,  K,  RHO).  Measurement  of  the  actual  boundary-layer  properties  should 
improve  the  accuracy  of  the  calculated  results. 
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The  accuracy  of  the  modified  inertia  parameter  is  quoted  as  +  5  percent  in  Ref.  6. 
Accuracy  may  be  improved  by  using  a  more  rigorous  water  droplet  trajectory  analysis.  This 
requires  an  accurate  description  of  the  flow  field  about  the  leading  edge  which  can  be 
obtained  from  experimental  measurements  or  estimated  using  computational  fluid 
dynamics.  If  the  leading  edge  being  evaluated  is  in  close  proximity  to  other  aerodynamic 
surfaces  (e.g.,  the  airframe  forebody,  control  surfaces,  etc.),  the  effects  of  the  changes  in 
flow  field  from  ideal  free-stream  conditions  should  be  accounted  for  in  the  model  inputs. 

Icing  tests  of  airframe  systems  have  become  more  common  in  recent  years  in  wind- 
tunnel-type  ground  test  facilities  because  of  improved  capabilities  for  producing, 
controlling,  and  measuring  icing  environmental  conditions.  However,  because  of  recent 
escalations  in  energy  costs,  wind  tunnel  test  time  has  become  increasingly  expensive.  Test 
time  can  be  minimized  by  predicting  the  effects  of  LWC  on  the  test  article  surface 
temperature  for  a  given  free-stream  static  temperature  or  vice  versa.  In  addition  measured 
data  can  be  input  to  the  model  and  the  results  used  to  analyze  and  evaluate  the  measured 
parameters.  Examples  of  model  applications  are  presented  in  the  remainder  of  this  section 
for  an  airfoil  leading-edge  thickness  of  0.045  ft.  The  validated  model  was  used  for  the 
calculations. 

5.1  EFFECTS  OF  VARYING  LIQUID  WATER  CONTENT 

5.1.1  Surface  Recovery  Temperature 

One  application  of  the  model  is  to  determine  the  average  equilibrium  surface 
temperature  (or  recovery  temperature)  of  the  leading  edge  for  a  given  set  of  environmental 
conditions.  As  an  example,  the  surface  temperature  was  calculated  as  a  function  of  LWC  for 
two  flight  Mach  numbers  at  3,000-ft  altitude.  Setting  TO  =  30°F  and  QK  and  PA  =  0.0,  the 
model  was  iterated  on  TS  until  QNET  =  0.0  ( |QNET|  <  0.01  w/in.2,  TS  within  ±0.05°F). 
The  results  for  flight  Mach  numbers  0.4  and  0.8  are  plotted  in  Fig.  24a.  At  0.4  Mach 
number,  the  average  surface  temperature  decreased  from  43.5°F,  the  dry  recovery 
temperature,  to  40.5°F  as  LWC  increased  from  0.0  to  about  0.013  gm/m3  or  a  rate  of 
change  of  230°F/gm/m3.  Further  increase  in  LWC  up  to  3.0  gm/m3  reduced  the  surface 
temperature  to  36.5 °F,  only  an  additional  4.0°F.  Over  the  full  range  of  LWC  change  from 
0.0  to  3.0  gm/m3,  the  total  change  in  surface  temperature  was  only  7.0°F.  Therefore,  at  0.4 
Mach  number,  TS  is  relatively  insensitive  to  LWC.  At  0.8  Mach  number,  however,  the  initial 
rate  of  change  was  higher  (304°F/gm/m3),  but  the  total  change  for  the  same  change  in  LWC 
was  28.4°F.  As  noted  previously,  the  inflection  points  occur  at  the  transition  of  MPR  from  a 
value  of  unity  to  a  value  less  than  unity  or  where  runback  initiates  as  LWC  is  increased. 
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A  typical  range  of  test  LWC’s  is  0.1  to  1.5  gm/m3.  For  this  range,  the  change  in  surface 
temperatures  at  flight  Mach  numbers  0.4  and  0.8  were  2.3  and  8.5°F,  respectively  (Fig.  24a). 
Since  this  15:1  change  in  LWC  produced  relatively  small  changes  in  surface  temperature, 
LWC  settings  in  this  range  are  not  critical  for  incipient  icing  tests. 

The  effects  of  flight  Mach  number  and  altitude  on  surface  temperature  were  investigated 
at  a  constant  LWC  =  1.5  gm/m3.  At  3,000-ft  altitude,  the  surface  temperatures  were  37.7 
and  57.3° F  for  0.4  and  0.8  Mach  numbers,  respectively,  or  a  difference  of  19.6°F  (Fig.  24a). 
Referring  to  Fig.  24b  for  0.4  Mach  number,  the  calculated  surface  temperatures  at  3,000-and 
1 5 ,000-ft  altitudes  (LWC  =  1.5  gm/m3)  were  37.8  and  36.9°F,  respectively,  a  difference  of 
only  about  1 .0°F.  Therefore,  surface  temperature  is  affected  more  by  flight  Mach  number 
than  by  altitude. 

5.1.2  Free-Stream  Static  Temperature 

Another  application  of  the  model  is  to  evaluate  the  relative  change  in  free-stream  static 
temperature,  TO,  with  LWC.  The  procedure  used  was  to  set  TS  =  34°F,  QK  and  PA  =  0.0, 
and  to  iterate  the  model  on  TO  until  QNET  «  0.0  (absolute  value  of  QNET  :£  0.01  w/in.2, 
as  before)  while  increasing  LWC  from  0.0  to  3.0  gm/m3  at  3,000-ft  altitude.  The  effect  of 
changing  flight  Mach  number  from  0.4  to  0.8  was  also  calculated,  and  the  results  are 
presented  in  Fig.  25a.  The  total  change  in  TO  for  an  increase  in  LWC  from  0.0  to  3.0  gm/m3 
was  6.6  and  25.7°F,  respectively,  for  flight  Mach  numbers  of  0.4  and  0.8. 

The  effect  of  changing  altitude  is  shown  similarly  in  Fig.  25b.  Again,  the  effect  of 
altitude  is  relatively  small  for  the  0.0  to  3.0  gm/m3  change  in  LWC.  The  changes  were  6.5 
and  7.2°F  at  altitudes  of  3,000  and  15,000  ft,  respectively. 

5.2  ICING  ONSET  PREDICTIONS 

Data  generated  by  iteration  of  the  model  on  TO  can  be  used  to  provide  predictions  of 
icing  onset  for  a  32°F  surface  temperature.  A  primary  test  objective  of  one  of  the  two  test 
configurations  discussed  previously  was  to  determine  the  inlet  leading-edge  icing  onset 
temperatures.  The  thermal  model  was  used  to  predict  icing  onset  for  a  range  of  flight 
conditions  prior  to  acquisition  of  the  test  data  used  for  model  validation.  The  free-stream 
and  local  flow  properties  were  assumed  identical,  as  in  the  case  of  an  isolated  airfoil,  one 
where  the  flow  field  is  unaffected  by  proximate  bodies.  The  pretest  predictions  for  a  32°F 
surface  temperature  are  shown  in  Fig.  26  for  three  altitudes  and  liquid  water  contents. 
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A  comparison  of  the  predictions  and  test  results  is  shown  in  Fig.  27.  The  upper  and  lower 
boundaries  identified  in  Fig.  26  are  replotted  in  Fig.  27  as  dashed  lines.  The  test  results  are 
represented  by  the  cross-hatched  band.  The  isolated  airfoil  predictions  agreed  with  the  test 
results  within  1  to  4°F  only  at  a  flight  Mach  number  of  0.6.  At  Mach  number  0.5,  icing 
occured  at  free-stream  static  temperatures  about  5°F  higher  than  predicted.  Above  Mach 
number  0.6,  the  predicted  temperatures  were  higher  than  the  test  results,  up  to  about  14°F  at 
Mach  number  0.8.  The  effect  of  using  the  0.3  evaporative  correlation  coefficient  (the 
validated  model)  is  shown  by  the  solid-line  boundaries.  This  improved  the  agreement  at  the 
higher  Mach  numbers  but  degraded  the  agreement  at  the  lower  flight  speeds.  A  more 
extensive  survey  of  the  local  flow  conditions  should  improve  the  results,  but  further  model 
development  is  also  needed  to  improve  the  icing  onset  predictions. 

5.3  ANTI-ICING  HEATER  SIZING 

The  model  can  also  be  used  for  determing  the  leading-edge  heater  capacity.  Two  such 
applications  are  discussed  below. 

5.3.1  Reducing  Icing  Onset  Temperature 

Referring  to  Fig.  26,  the  family  of  curves  can  be  shifted  to  the  left  by  applying  internally 
generated  heat  to  the  leading  edge  using,  for  example,  a  heating  device,  PA,  as  shown 
schematically  in  Fig.  1 .  To  illustrate  the  effect  of  heat  addition,  the  iterative  solution  for  the 
icing  onset  boundary  for  TS  =  32°F,  3,000-ft  altitude,  and  LWC  =  0.5  gm/m3  is  shown  in 
Fig.  28  (QK  =  0.0).  The  curve  was  shifted  to  the  left  by  incrementally  setting  PA  at  5  and  10 
w/in.2,  decreasing  the  free-stream  static  temperature  at  which  icing  onset  would  occur  at  0.5 
flight  Mach  number  by  18.6  and  38.5°F,  respectively.  This  technique  can  be  used  along  with 
the  flight  system  mission  objective  to  improve  the  margins  for  ice-free  operation,  i.e.,  to 
enlarge  the  ice-free  flight  envelope. 

5.3.2  Dry  Anti-Icing 

For  some  flight  systems,  it  may  be  desired  to  size  a  heater  for  complete  evaporative  anti- 
icing,  referred  to  by  Sogin  (Ref.  9)  as  dry  anti-icing.  To  satisfy  this  requirement,  TS  must  be 
increased  until  MPR  =  1.0,  i.e.,  until  all  of  the  intercepted  water  is  evaporated  (zero 
runback).  Application  of  the  validated  model  for  dry  anti-icing  is  illustrated  in  Fig.  29  for 
three  liquid  water  contents  at  the  0.4  Mach  number,  3000-ft  altitude,  20°F  TO  flight 
condition.  The  power  density  required  is  8.9,  32.7,  and  58.6  w/in.2,  respectively,  for  liquid 
w'ater  contents  of  0. 1 , 0.5,  and  1 .0  gm/m3.  The  dry  anti-icing  power  density  was  obtained  by 
increasing  TS  until  MPR  =  1.0.  The  corresponding  surface  temperatures  required  for  dry 
anti-icing  are  67.2,  1 14.0,  and  138. 3°F. 
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5.4  TEST  DATA  ANALYSIS 

The  thermal  model  can  also  be  used  to  analyze  and  evaluate  the  leading-edge  thermal 
balance  data  obtained  from  wind  tunnel  test  and/or  flight  test.  Careful  attention  should  be 
given  to  the  uncertainty  of  the  measured  inputs  XMO,  PO,  LWC,  DD,  TS,  D,  A,  QK,  and 
PA.  Some  insight  on  this  subject  can  be  obtained  from  the  sensitivity  study  presented  in 
Section  4.1  of  this  report.  LWC  and  DD  require  calibration  of  the  wind  tunnel  water  spray 
system  or  measurement  during  the  flight  test.  The  leading-edge  thickness,  D,  and  surface 
area,  A,  require  the  measurement  of  physical  dimensions  which  are  more  difficult  for  small¬ 
thickness  airfoils.  If  the  leading-edge  contour  is  not  strictly  a  half-cylinder,  then  the  surface 
area  should  be  measured  and  converted  to  an  equivalent  cylindrical  diameter  (using 
measurement  of  the  span),  or  the  more  appropriate  E  versus  KO  relationship  for  the 
particular  leading-edge  shape  should  be  used. 

For  test  data  analysis,  the  model  should  be  calibrated  or  verified.  The  ideal  procedure 
would  be  to  check  each  term  of  the  heat  balance  equation  individually.  Since  the  terms  are 
interactive,  however,  the  terms  can  only  be  checked  in  groups.  This  can  be  done 
progressively  in  three  steps.  By  inputting  the  measured  parameters  for  a  dry-air  calibration 
test  (Step  1),  the  terms  QC  and  QF  can  be  evaluated  quantitatively  as  a  group,  Similarly, 
during  a  wet-air  calibration  test  (Step  2),  the  effect  of  the  terms  which  are  functions  of  LWC 
(QE,  QS,  and  QKE)  plus  QC  and  QF,  all  five  terms,  can  be  investigated  as  a  group.  And 
finally,  with  the  heater  on  (Step  3),  the  computed  results  for  the  complete  equation  should 
result  in  thermal  balance. 

Steps  1  and  2  require  that  no  (or  negligible)  thermal  gradients  exist  across  the  structural 
interface.  This  must  be  confirmed  by  measurement.  For  Step  3,  measurement  of  the  gradient 
is  also  required  since  a  temperature  difference  across  the  interface  will  probably  exist  with 
the  heater  on. 

Step  3  requires  an  input  for  QK.  Quantitative  calculation  of  QK  is  usually  difficult 
because  the  structural  complexity  of  the  interface  configuration  complicates  the 
determination  of  the  effective  conduction  heat-transfer  coefficient  and  because  materials  are 
often  used  for  which  the  coefficient  is  unknown.  However,  QK  can  be  estimated  using  the 
following  procedure.  First  (Step  A),  measure  the  temperature  gradient  (AT)  across  the 
structural  interface  while  flowing  dry  air.  With  QNET  set  equal  to  zero  (also  QE,  QS,  and 
QKE  =  0.0  for  dry  air)  and  the  heater  off,  then  QK  equals  0.0  (provided,  of  course,  that 
AT  *=  0.0)  and  QC  equals  QF.  Secondly  (Step  B),  repeat  Step  A  with  the  heater  on  at 
constant  power,  again  measuring  AT.  Assuming  that  QF  still  equals  QC,  then  the  following 
cancellation  applies: 
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QK  =  jQT  +  PA  -J$e'=  (K/X)  AT 

Therefore,  the  conduction  heat-transfer  coefficient,  K,  equals  (PA)X/AT.  Step  B  should  be 
repeated  for  each  anticipated  heater  power  setting. 

Steps  3  and  B  also  require  an  input  for  PA.  The  most  convenient  arrangement  for 
determining  the  power  density  is  the  use  of  an  electrical  heater  for  which  the  electrical  power 
consumption  is  measured  directly.  Heater  end  fittings  or  other  extraneous  sources  of  heat 
loss  require  close  scrutiny. 

Model  verification  Steps  1  through  3  should  be  accomplished  for  each  planned  flight 
condition. 

Once  the  above  steps  have  been  completed,  one  can  evaluate  the  thermal  characteristics 
of  the  leading-edge  system  by  inputting  the  measured  parameters  and  plotting  QG  versus 
QL.  Ideally,  the  data  should  fall  on  a  line  connecting  all  points  where  QG  -  QL  or  QNET 
=  0.0.  In  reality,  QNET  *  0.0.  Some  assessment  of  the  TS  data  quality  can  be  made  by 
omitting  TS  as  an  input  (since  the  model  is  most  sensitive  to  TS),  iterating  the  model  on  TS 
until  QNET  -  0.0,  and  evaluating  the  differences  between  TS  calculated  and  TS  measured. 

6.0  SUMMARY  OF  RESULTS 

1.  A  relatively  simple  steady-state  thermal  model  has  been  developed  to  assist  in  predicting 
incipient  icing  on  wings,  fins,  and  inlet  leading  edges  of  small  subsonic  missiles  and 
drones. 

2.  The  model  is  primarily  a  heat  balance  of  the  component  leading  edge  and  treats  the 
effects  of  free-stream  Mach  number,  altitude,  cloud  liquid  water  content,  and  vehicle 
anti-ice  heaters. 

3.  The  model  has  proved  to  be  useful  in  predicting  the  effects  of  liquid  water  content  on 
vehicle  surface  temperature,  hence  the  onset  of  icing,  during  free-jet  wind  tunnel  tests. 
This  capability  has  reduced  the  required  test  time,  thus  reducing  energy  costs  and  total 
test  costs. 

4.  The  model  is  useful  for  related  applications,  e.g., 

a.  Determining  the  free-stream  static  temperature  for  a  specified  surface  temperature 
with  varying  liquid  water  content 
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b.  Sizing  anti-icing  heaters 

c.  Analyzing  wind  tunnel  or  flight  test  incipient  icing  test  data. 

5.  Analytical  application  of  the  steady-state  model  has  been  limited  to  the  following 
conditions: 

a.  Leading-edge  thickness,  0.22  to  0.54  in. 

b.  Volume  median  droplet  diameters,  10  to  30  n- 

c.  Altitudes  up  to  15,000  ft 

d.  Mach  numbers  up  to  0.8 

e.  Free-stream  static  temperatures,  -25  to  32°F 

f.  Leading-edge  surface  temperatures,  32  to  130°F 

g.  Liquid  water  contents  up  to  3.0  gm/m3. 

6.  Using  the  available  AEDC  test  data,  correlations  of  the  calculated  and  measured 
leading-edge  surface  temperatures  were  used  to  validate  the  thermal  model.  The  average 
model  bias  is  about  +4°F  at  a  surface  temperature  of  32°F. 
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Steadv-State  Heat  Balance  Aquation 


QNET 


QC  QE  QS 

=  HC  ( TS  -  TO)  +  (EC)*  MPR  ( M  /  A )  L  +  (M/A)  CW  (TS  -  TO)  +  QK 

Convective  Evaporative  Sensible  Conduction 

QF  QKE 

-  HC ( R )  V2/2GJ(CA)  -  (M/A)  V2/2GJ  -  PA 

Air  Water  Heater 

Kinetic  Energy 


Heat 

T.nss 

Terms 


Heat 

Cain 

Terms 


*EC ,  Evaporative  Correlation  Coefficient  (Added  During  Model  Validation) 

Figure  1.  Typical  leading-edge  system  and  heat  balance  equation. 
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NOTES 

J  refers  to  Joukowski  0015  airfoil. 

Subscript  c  means  cambered. 

*Curve  for  cylinder  used  to  determine 
water  catch  efficiency,  E. 

This  figure  taken  from  Ref.  6. 
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MODIFIED  INERTIA  PARAMETER,  K0 

Figure  2.  Water  catch  efficiencies. 
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Heat  Gain  (QG),  w/in. 


m 


T  0.3  w/m 


QNET  =  0.0  (Desired 
Thermal  Balance 
or  Line  of  Ideal 
Agreement) — ..  ^ 


Fi rst-Order 
Curve  Fit 


f  C£\£&  ^2.0 
£  w /  in  •  2 


w 


SD  =  +1.23  w/in. 


O  Heater  Off 
A  Heater  On 


Heat  Loss  (OT,),  w/in.2 

Figure  3.  Thermal  balance,  dry  air  (Configuration  A). 
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Surface  Temperature  Difference 
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Measured  Surface  Temperature  (TSM),  °F 

Figure  S.  Difference  between  measured  and  calculated  leading-edge  surface 
temperatures,  dry  air  (Configuration  A). 
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Measured  Surface  Temperature  (TSM),  °F 

Figure  7.  Calculated  versus  measured  surface  temperature,  wet  air  (Configuration  A) 


Heat  Loss  (QL),  w/in. 

Figure  8.  Thermal  balance,  model  adjusted,  wet  air  (Configuration  A) 


Curve  Fit 


Measured  Surf  a  tie  Temperature  (TSM),  °F 

Figure  9.  Calculated  versus  measured  surface  temperature, 
model  adjusted,  wet  air  (Configuration  A), 


Heat  Gain  (QG) ,  w/in. 
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Unbalance  in  QNET  (or  QK),  w/in. 
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Figure  13.  Thermal  balance,  dry  air,  QK  =  f  (DT)  (Configuration  B). 
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Calculated  Surface  Temperature  (TSC),  °F 
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Calculated  Surface  Temperature  (TSC),  °F 
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Figure  17.  Calculated  versus  measured  surface  temperature,  wet  air  (Configuration  B). 
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Heat  Loss  (QL) ,  w/in.2 


Figure  18.  Thermal  balance,  model  adjusted,  wet  air  (Configuration  B). 
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Figure  19*  Calculated  versus  measured  surface  temperature,  model  adjusted,  wet  air 
(Configuration  B). 
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Convective  Ht 

Convective  Ht  Fraction  of  Water  Water  Catch  Modified  I  nerti3  Transfer  Evaporation 

Transfer  (QC)  Evaporated  (MPR)  Rate  1M(A}  Parameter  (KO)  Coefficient  (HC)  Potential  TEP) 
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A 

Leading-Edge  Area 

S 

Airfoil  Span 

D 

Leading-Edge  Effective  Diameter 

PO 

free-Stream  Static  Pressure 

DD 

Droplet  Diameter 

TO 

Free-Stream  Static  Temperature 

LWC 

Liquid  Water  Content 

TS 

Leading-Edge  Surface  Temperature 

PA 

Heater  Power  Density 

XMD 

Flight  Mach  Number 

Figure  21.  Sensitivity  of  output  parameters  for  a  one-percent 
increase  in  input  parameters. 
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Change  in  QNET ,  w/in. 
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Figure  22.  Sensitivity  of  QNET  for  a  one-percent 
increase  in  input  parameters. 
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Figure  23.  Typical  heat  loss  audit  as  a  function  of  liquid  water  content. 
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TO  =  30°F 
QK,  PA  -  0.0 
D  -  0.045  ft 
DD  =  20.0  y 


Fraction  of  Water 
Evaporated  (MPR) 

-  <  1.0 

_  1.0 


0  0.02  0.04  0.06  0.08  0.08  1.0  2.0  3.0 

Liquid  Water  Content  (LWC),  gm/m3 

a.  3,000-ft  altitude 

Figure  24.  Effect  of  liquid  water  content  on  leading-edge  surface  temperature. 
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Surface  Temperature 


Fraction  of  Water 

TO  =  30°F  Evaporated  (MPR) 

QK,  PA  -  0.0 

D  -  0.045  ft  - <  1.0 


b.  0.4  Flight  Mach  No. 
Figure  24.  Concluded. 
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Free-Stream  Static  Temperature 
(TO),  of 


0  0.01  0.02  0.03  0.03  1.0  2.0  3.0 

Liquid  Water  Content  (LWC),  gm/m^ 

a.  3,000-ft  altitude 

Figure  25.  Effect  of  liquid  water  content  on  free-stream  static  temperature. 
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Free-St ream  Static 
Temperature  (TO),  ° 


TS  =  S4°F 
QK,  PA  -  0.0 
D  =  0.045  ft 
DD  =  20.0  u 


Fraction  of  Water 
E  vapor  a  ted  (  MPR ) 

-  <  1.0 

-  1.0 


b.  0.4  Flight  Mach  No. 
Figure  25.  Concluded. 
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Flight  Mach  Number  (XMO) 


Figure  26.  Icing  onset  predictions,  isolated  airfoil  (Configuration  B). 
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Figure  27.  Comparison  of  icing  onset  predictions  and  test  results  (Configuration  B). 
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Figure  28.  Effect  of  heater  on  reducing  icing  onset  temperature. 
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Table  1.  Uncertainties  of  Calculated  Data 

3,000-ft  Altitude,  0.55  Mach  No.,  0.8-gm/m3  Liquid  Water  Content 

Uncertainty,  Heat  Uncertainty,  Heat 

Data  Gain  (QG),  percent  Loss  (QL),  percent 

Configuration  A  4.72  9.05 

Configuration  B  7.60  5.25 
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1 -percent 

Parameter  Base  Increase  In 

Units  Values  Input 

Parameter 


m  20.00  20.20 

gm/m3  0.80  0.808 

w/in.2  5.00  5.05 

psia  13.171  13,30 


°F  34.00  34.34 

—  0.55  0.5555 

ft2  0.16  0.1616 

ft  0.045  0.04545 

ft  2.30  2.323 
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NOMENCLATURE 

A  Leading-edge  surface  area,  ft2 

AO,  A1  First-order  curve-fit  coefficients 
ALT  Geopotential  altitude,  ft 

CA  Specific  heat  at  constant  pressure  for  air,  Btu/lbm-°F 
CW  Specific  heat  of  water,  Btu/lbrn-^T 

D  Leading-edge  thickness  or  effective  diameter,  in.,  ft 

DD  Volume  median  water  droplet  diameter,  n  (microns) 

DF  Diffusion  coefficient,  ft2/hr 

DT  Interface  temperature  gradient,  °F 

E  Water  catch  efficiency 

EC  Evaporative  correlation  coefficient 

EP  Evaporation  potential,  lbm/hr-ft2 

G  Constant,  32,174  lbm-ft/lbf-sec2 

HC  Calculated  convective  heat-transfer  coefficient,  Btu/hr-ft2-°F 

HT  Heat 

1  Dimensionless  parameter 

J  Mechanical  equivalent  of  heat,  778  ft-lbf/Btu 

K  1,000;  also  conduction  heat-transfer  coefficient,  Btu/hr-ft-°F 

KM  Mass  transfer  function,  lbm-°F/Btu 
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KO  Modified  inertia  parameter 

L  Latent  heat  of  vaporization  for  water,  Btu/lbm 

LWC  Liquid  water  content,  gm/m3 

M  Water  catch  rate,  lbm/hr 

M/A  Water  catch  rate  per  unit  area,  lbm/hr-ft2 

MPR  Fraction  of  water  evaporated 

MU  Air  viscosity,  Ibf-sec/ft2 

PA  Heater  power  density,  Btu/hr-ft2,  w/in.2 

PO  Free-stream  static  pressure  (ambient),  psia 

PR  Prandtl  number 

PTO  Free-stream  total  pressure,  psia 

PTOW  Water  vapor  pressure  at  TO,  psia 

PTSW  Water  vapor  pressure  at  TS,  psia 

QC  Convective  heat  transfer,  Btu/hr-ft2,  w/in.2 

QE  Evaporative  heat  transfer,  Btu/hr-ft2,  w/in.2 

QF  Air  kinetic  energy,  Btu/hr-ft2,  w/in.2 

QG  Heat  gain,  sum  of  QF,  QKE,  and  PA,  Btu/hr-ft2,  w/in.2 

QK  Conductive  heat  transfer,  Btu/hr-ft2,  w/in.2 

QKE  Water  kinetic  energy,  Btu/hr-ft2,  w/in.2 

QL  Heat  loss,  sum  of  QC,  QE,  QS,  QK,  Btu/hr-ft2,  w/in.2 
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QNET  Net  heat  transfer,  QG  -  QL,  BTU/hr-ft2,  w/in.2 
QS  Sensible  heat  transfer,  Btu/hr-ft2,  w/in.2 

R  Boundary-layer  recovery  factor 

RHO  Air  density,  lbm/ft3 

S  Airfoil  span  (inlet  perimeter),  ft 

SC  Schmidt  number 

SD  Estimate  of  the  one-sigma  standard  deviation 
Tl,  2, 

etc.  Temperature,  °F 

TF  Reference  temperature,  °R 

TFK  Reference  temperature,  °K 

TO  Free-stream  static  temperature  (ambient),  °F 
TTO  Free-stream  total  temperature,  °F,  °R 

TS  Leading-edge  surface  temperature,  °F,  °R 

TSC  Calculated  surface  temperature,  °F 

TSM  Measured  surface  temperature,  °F 

V  Flight  velocity,  ft/sec 

VK  Flight  velocity,  knots 

X  Average  thickness  of  structural  interface,  ft 

XMO  Flight  Mach  number 
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